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1 Introduction 



Measure-valued branching processes or superpro cesses constitute a rich class of infinite 
dimensional processes currently under rapid development. Such processes arose in appli- 
cations as high density limits of branching particle systems; see e.g. Dawson (1992, 1993), 
Dynkin (1993, 1994), Watanabe (1968). The development of this subject has been stimu- 
lated from different subjects including branching processes, interacting particle systems, 
stochastic partial differential equations and non-linear partial differential equations. The 
study of superprocesses has also led to a better understanding of some results in those 
subjects. In the literature, several different types of supcrproccss have been introduced 
and studied. In particular, Dawson and Hochbcrg (1991), Dawson et al (1990) and Wu 
(1994) studied multilevel branching superprocesses, Gorostiza and Lopez- Mimbela (1990), 
Gorostiza and Roelly (1991), Gorostiza et al (1992) and Li (1992b) studied multitype 
superprocesses, Dynkin (1993, 1994) and Li (1992a, 1993) studied non-local branching 
superprocesses, Gorostiza (1994) studied mass-structured superprocesses. Hong and Li 
(1999) and Li (2002) studied superprocess-controUed immigration processes, and Bose 
and Kaj (2000) studied age-reproduction-structured superprocesses. Those models arise 
in different circumstances of application and are of their own theoretical interests. 

In this paper, we provide a unified treatment of the above models. We first give a 
new formulation of the non-local branching superprocess as the high density limit of some 
specific branching particle systems. Then we derive from this superprocess the multitype, 
the mass-structured, the multilevel and the age-reproduction-structured superprocesses 
and superprocess-controUed immigration processes. Another related model, the so-called 
rebirth superprocesses, is also introduced to explain the non-local branching mechanism. 
This unified treatment simplifies considerably the proof of existence of the old classes of 
superprocesses and also gives rise to some new ones. We think that this treatment may 
give some useful perspectives for those models. The unification is done by considering 
an enriched underlying state space E x I instead of E. In this way, the mutation in 
types of the offspring can be modeled by jumps in the /-coordinates so that the multitype 
superprocess can be derived. The superprocess-controUed immigration process is actually 
a special form of the multitype superprocess. To get the mass-structured superprocess 
we let / = (0, oo), which represents the mass or size of the infinitesimal particles. For 
the age-reproduction-structured superprocess, we take / = [0, oo) x N, where N is the 
set of non-negative integers, to keep the information on ages and numbers of offspring 
of the particles. In this model we have of course that any particle has non-decreasing 
[0, oo) X N-coordinates starting from (0, 0) at its birth time. To get a two level super- 
process, we simply assume that / = M{S)° is the space of nontrivial finite measures on 
another space S, and the M(5') "-coordinate of the underlying process is a superprocess 
itself. For two-level branching systems, what has been done so far for the second level 
branching is local, that is, when a superparticle branches, superoffspring are produced 
as exact copies of their parent. Since the superparticles have an internal dynamics and 
evolve as branching systems themselves, it is desirable to have the possibility that the 
superoffspring have internal structures different from those of their parents, which re- 
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quires a non-local branching mechanism. Models of this type have potential applications 
in genetics, population dynamics and other complex multilevel systems; see e.g. Dawson 
(2000) and Jagers (1995). 

Notation and basic setting: Suppose that is a Lusin topological space, i.e., a homc- 
omorph of a Borel subset of some compact metric space, with Borel o"-algebra lS§{E). Let 
M{E) denote the space of finite Borel measures on E topologized by the weak conver- 
gence topology, so it is also a Lusin topological space. Let N{E) be the subspace of M{E) 
consisting of integer-valued measures on E and let M{E)° — M{E) \ {0}, where denotes 
the null measure on E. The unit mass concentrated at a point x & E is denoted by 5x- 
Let 

B{E) — { bounded =^(£J)-measurable functions on E }, 
C{E) = { f: f e B{E) is continuous }, 
Ba{E) = { /: / e B{E) and ||/|| < a }, 

where a > and "|| • ||" denotes the supremum norm. The subsets of positive members 
of the function spaces are denoted by the superscript e.g., B~^[E), C^{E). For 

/ e B{E) and fi e M(E), we write /x(/) for /d/^. 

2 Non-local branching particle systems 

Non-local branching particle systems have been considered by many authors. We here 
adapt the model of Dynkin (1993). Let ^ = (i7, ^t, Px) be a right continuous strong 

Markov process with state space E and transition semigroup {Pt)t>G- Let 7 £ B'^(E) and 
let F{x, di') be a Markov kernel from E to N{E) such that 



A branching particle system with parameters (^, 7, F) is described by the following prop- 
erties: 

(2. A) The particles in E move randomly according to the law given by the transition 
probabilities of ^. 

(2.B) For a particle which is alive at time r and follows the path {Ct)t>r, the conditional 
probability of survival during the time interval [r,t] is p{r,t) := exp{— J^. ^{^s)ds}. 

(2.C) When a particle dies at a point x E E, it gives birth to a random number of 
offspring in E according to the probability kernel F^Xjdv). The offspring then start to 
move from their locations. (Thus the name "non-local branching" is used.) 

In the model, it is assumed that the migrations, the lifetimes and the branchings of 
the particles are independent of each other. Let Xf(B) denote the number of particles in 
B G ^{E) that are alive at time t > and assume Xq{E) < 00. Then {Xf : t > 0} is 
a Markov process with state space N{E). For a G N{E), let denote the conditional 
law of {Xt:t> 0} given Xq = a. For / G B+{E), put 




(2.1) 



Ut{x) = Ut{x,f) 



logQs^eM-Mf)}- 



(2.2) 
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The independence hypotheses imply that 

Q<xexp{-Xt(/)} = exp{-a{ut)}. 
Moreover, we have the following fundamental equation 



(2.3) 



e-"'(^) = P.{p(0,t)e-^«')} + pJ / fp(0,s)7(e 



N{E) 



F(e„dz/)]ds}. (2.4) 



This equation is obtained by thinking that if a particle starts moving from point x at 
time 0, it follows a path of ^ and does not branch before time t, or it first splits at time 
s G {0,t]. By a standard argument one sees that equation ()2.4j) is equivalent to 



-ut{x) 



P^e-^(«') - P,.| ^ 7(e.)e-"'-(«=)d.s} 



(2.5) 



-U{ut-s) 



N{E) 



F(e.,dz/)]ds}; 



see e.g. Dawson (1992, 1993) and Dynkin (1993, 1994). It is sometimes more convenient 
to denote 



vt{x) = vt{x, /) = 1 - exp{-Mt(x)} 
and rewrite ()2.5|1 into the form 



VtiX 



(2.6) 



(2.7) 



{1 - e--^«')} - P,.{ ^ 7(es)^^t-s(e.)ds} 
+ PJ f U{Q [ (l-e-^('^'-»))F(e.,di/)ld4. 

Jo L Jn(E) ^ 



3 Non-local branching superprocesses 

In this section, we prove a limit theorem for a sequence of non-local branching particle 
systems. Although the particle systems considered here are very specific, they lead to 
the same class of non-local branching superprocesses constructed in Dynkin (1993, 1994) 
and Li (1992a) with a slightly different formulation. We shall give some details of the 
derivation to clarify the meaning of the parameters, which is needed in understanding the 
connections of non-local branching with other related models. 

Let {Xt{k) : t > 0}, = 1, 2, . . . be a sequence of branching particle systems with 
parameters (^,7fc,-Ffc). Then for each k, 

{Xf ) := k-^Xt{k) :t>0} (3.1) 

defines a Markov process in Mk{E) := {k'^a : a e N{E)}. For a E Mk{E), let Q^^^ 
denote the conditional law of {xj:''^ '■ t > 0} given X^^^ = a. By (Q we have 



Q«exp{-X«(/)}=exp{-a(A;nf)) 



(3.2) 
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where u[''\x) is determined by 



vl''\x) = k[l — exp{—u\'''{x)}] 



and 



P. {k{l - e-^«*)/'^)} - P,.{ ^ 7.fe)tSfe)ds~ 

+ PJ r [A;7fc(e.) / (1 -e-^("*-))Ffc&,d^/)ld4• 
'- Jo '- JAffB) -I 



(3.3) 



(3.4) 



For fi G M{E), let 0"^^ be a Poisson random measure on E with intensity kfi, and let Q^'^^ 



denote the conditional law of {xj:'''^ '■ t > 0} given 

^0 = ^"'^fcM- From dH we get 



(m) 



Q|3exp{-Xf (/)} = exp{-M^^f)}. 



(3.5) 



It is natural to treat separately the offspring that start their motion from the death 
sites of their parents. Suppose that G B^{E x [0, 1]) and, for each x E E, 

oo 

gk{x,z) = Y,pt\^)^\ 2;g[0,1], 

1=0 

is a probability generating function with snp^(zE{d/dz)gk{x,l^) < oo. Let and Pk G 
B^{E) and assume 7fc(a;) := ak{x) + I3k{x) is strictly positive. Let F^^^x^du) be another 
probability kernel from E to N{E) satisfying ()2.H) . We may replace Ffc(a;,di/) by 



7k{x) 



ak(x^ 



(3.6) 



i=0 



where Flf\x,di') denotes the unit mass concentrated at i6x- Intuitively, as a particle 
splits at a; G -E, the branching is of local type with probability ak{x)/'yk{x) and is of 
non-local type with probability l3k{x)/'-fk{x)- If it chooses the local branching type, the 
distribution of the offspring number is {p[''\x)}. The non-local branching at x G -E is 
described by the kernel F^''\x,di'). Now ()3.4|1 turns into 



(k) 
vl [X 



) = Px.{A;(l 



)} - P.{ l\kiQvt\iQds] 



+ Px 



I kakiQ[l - gkiCs,e~ 



(3.7) 



kPkiis 



f [l-e-^("*-)]F(^)(e„dz/)ld4, 
Jn{e) ^ ^ 



or equivalently 

vl'\x) + f P.[0,(e., vf\{is)) + Mis, vfX)]ds = P.fc[l - e 
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_ p-/(?t)/fci 



(3.8) 



where 

(f>k{x,z) = kak{x)[gk{x, 1 - z/k) - (1 - z/k)] (3.9) 

and 

M^J)=Pk{x)[f{x)-Ck{x,f)], (3.10) 

where 

Ckix, f)= [ kil - exp{z/(log(l - f/k))})F^'\x, du). (3.11) 
Jn{e) 

Let Mq{E) denote the set of all Borel probability measures on E. Suppose that 
hk e B+{E X Mo{E) X [0, 1]) and, for each {x,tt) e Ex Mq{E), 

oo 

hk{x,7T,z) = ^ql''\x,7r)z\ z G [0,1], 

is a probability generating function with sup^^{d/dz)hk{x,n,l^) < oo. Suppose that 
G{x, dvr) is a probability kernel from E to Mq{E). We may consider a special form of the 
second term in ()3.6|) by letting 

„ oo 

F^^\x,diy)= [Vgf^(x,7r)(/7r)«(dz/)lG(x,d7r), (3.12) 

where Z7r(di^) denotes the image of vr under the map y ^ Sy from E to M{E) and (/vr)** 
denotes the i-fold convolution of In. Now we have 



Ckix, f) = k[l- hk{x, vr, 1 - 7T{f)/k)]G{x, dvr). (3.13) 

JMo{E) 

Intuitively, if a parent particle at x E E chooses non-local branching, it first selects 
an offspring-location-distribution tt{x, ■) G Mq{E) according to the probability kernel 
G{x, dvr), then gives birth to a random number of offspring according to the distribution 
{g-'^^(x, 7r(x, ■))}, and those offspring choose their locations in E independently of each 
other according to tt{x, ■). A similar non-local branching mechanism was considered in Li 
(1992a, 1993). 

In view of ()3.5|) and ()3.8p . it is natural to assume the sequences {(pk}, {Pk} and {(k} to 
converge if one hopes to obtain convergence of {xj:'^^ : t > 0} to some process {Xt : t > 0} 
as A; — ^ oo. 

Lemma 3.1 (i) Suppose that 

oo 

J2^qt\x,7r)<l (3.14) 

i=0 



and that Ck{x,f) — > C{x,f) uniformly on E x B^{E) for each a > 0, then C{x,f) has 
representation 

C(x, /) = Xix, f)+ [ (1 - e-^(^))r(a:, dz/), (3.15) 

Jm(e)° 

where X{x, dy) is a bounded kernel on E, and z/(l)r(x, du) is a bounded kernel from E 
to M{EY with 

X{x,l)+ I u{l)r{x,du) <1. (3.16) 

(a) A functional ({x, f) can be given by i3.15}) and (ltj.J6]) if and only if it has repre- 
sentation 

r r 1 

ax,f)= dix,7i)7rif)+ (l-e-"-(^))n(a:,7r,dM) G(x,d7r), (3.17) 

J Mo{E) L JQ ^ 

where d G B'^i^E x Mq[E)), un{x, n, du) is a bounded kernel from E x Mo{E) to (0, oo) 
and G{x, dvr) is a probability kernel from E to Mq{E) with 

POD 

d{x,7T)+ / un{x,'K,du) < 1. (3.18) 
Jo 

(Hi) To each function C('5 ■) given by i3.17\) and l\3.1^) there corresponds a sequence 
of the form H3.1!^) satisfying the requirement of (i). 

Proof, (i) Note that k{l - converges to / uniformly in B^{E). Then 

Ck{x,k{l-e-^/'))= [ k{l-exp{u{f)/k})F^''\x,du) 

J N(E) 

converges to C(x, /) uniformly on E x B^ {E) . It is known that a metric r can be introduced 
into E so that [E, r) becomes a compact metric space while the Borel a-algebra induced 
by r coincides with ^{E)\ see e.g. Parthasarathy (1967, p. 14). Now M{E) endowed with 
weak convergence topology is a locally compact metrizable space. Let M{E) = M{E) U 
{d} be the one-point-compactification of M{E). By ()3.14|) . {A;z/(1)F('')(x, d(A;"V)) : x G 
E,k > 1} viewed as a family of finite measures on M{E) is tight. Fix x & E and take 
{ui] C {n} such that /cjZ/(l)Ffc.(x, d(fcj~^z/)) converges to some finite measure G{x,dh') on 
M{E) as i — >• cxD. It follows that 

C(x,/)= / {l-e-'''~f^)u{l)-'G{x,du), 

J M{E)° 

first for / G C^{E, r) and then for all / G B^{E). Now ()3.15p follows by a simple change 
of the measure and ()3.16p follows from ()3.14j) . (ii) is immediate. To get (iii) we may set 

POO 

hk{x, TT, z) = 1 + d{x, Tx){z - 1) + k-^ / (e^"("-i) - l)n(x, vr, du). 

Jo 
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Observe that 



—^hk{x,7r,0)>0, i = l,2, 

and ()3.14p assures that hk{x,n,0) > 0. Thus for fixed (x, a) E E x Mq{E), hk{x,Tr,-) is 
a probabihty generating function. Then we define Ck{x,f) by (j3.13p so that Ck{.x, f) = 
C(x,/)for(x,/)Gi?x5+,(E). □ 

Lemma 3.2 (Li, 1992c) (i) Suppose that, for each I > 0, the sequence (j)k{x,z) is uni- 
formly Lipschitz in z on the set E x [0,/] and that (j)k{x,z) converges to some z) 
uniformly as k —>■ oo, then (j){x, z) has the representation 

POO 

(f){x,z) = b{x)z + c{x)z'^ + (e"^"- l + zM)m(x, dw), xeE,z>0, (3.19) 

^0 

where b G B{E), c G B^{E) and {u A v?)m{x, du) is a bounded kernel from E to (0, oo). 

(a) To each function 0(-, ■) given by \3.iy\) there corresponds a sequence of the form 
\3.y\) satisfying the requirement of (i). 

Based on Lemmas 13.11 and 13.21 the following result can be proved similarly as in 
Dawson (1992, 1993), Dynkin (1993, 1994) and Li (1992a, c). 

Lemma 3.3 If the conditions of Lemma \3.1\ (i) and Lemma \3.2i (i) are fulfilled and if 

f3k (3 & B~^{E) uniformly as k ^ oo, then for each a > both v^''\x, /) and kuf^\x, /) 
converge boundedly and uniformly on the set [0, a]x E x B^{E) of {t, x, f) to the unique 
bounded positive solution Vtf{x) to the evolution equation 

Vtfix) + £ { lj<P{y, Vt-sfm + ^(y, Vt.sf)]Ps{x, d?/)}ds = Ptf{x), t > 0, (3.20) 
where 

^(x, /) = Pix)[fix) - C(x, /)], xeEJe B+iE). (3.21) 
By Lemma (3.31 and Dawson (1993, p. 42), 

/ e-^(^)Qi(/i, du) = exp{-/i(rj)}, / G 5+(i?), (3.22) 
Jm{e) 

defines a transition semigroup {Qt)t>o on M{E). A Markov process {Xt : t > 0} with 
state space M{E) is called a non-local branching superprocess with parameters (^, 0, ip) 
if it has transition semigroup {Qt)t>o- Condition ()3.14|) means that the corresponding 
branching particle system has subcritical non-local branching. In terms of the limiting 
superprocess, this condition is expressed as ()3.16|) . which is of course a restriction of the 
class of ■) given by (|3.15p . However, since 0(x, z) + [3{x)z belongs to the class defined 
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by ()3.19|) . and since P G B^{E) is arbitrary, ()3.16|) does not put any restriction on the 
generality of 

f{x)) + /) = /(x)) + (3{x)f{x) - /3(x)C(x, /). 

Therefore, the class of non-local branching superprocesses given by ()3.20j) and ()3.22p 
coincides with those constructed Dynkin (1993, 1994) and Li (1992a), where the first 
term of il){x,-) was written into 0(x, ■). In principle, (j3.2Uj) and (j3.22j) give the most 
general non-local branching superprocesses constructed in the literature up to now. A 
more general class of non-local branching superprocesses were discussed in Dynkin et al 
(1994), but their existence has not been established. The next theorem follows similarly 
as in Li (1992a, c). 

Theorem 3.1 Let {xf ^ : t > 0} be the sequence of renormalized branching particle 
systems determined by \3.5}) and and let {Xt : t > 0} be the non-local branching 

superprocess with transition semigroup {Qt)t>o given by i3.2(Jp and i3.22p . Assume that 
the conditions of Lemma \3.1\ (i) and Lemma \3.2l are fulfilled. Then for every fi G M{E), 
< ti < . . . < tn and a >0, as k ^ CO, 

n n 

1=1 i=i 

uniformly on fi, . . . Jn ^ B+{E). 

Naturally, we may regard {X^'^^ : t > 0} as a process in the space M{E). Then the 
above theorem shows that the finite dimensional distributions of {X^'^^ '■ t > 0} under 
converge as ^ oo to those of {Xt : t > 0} under Q^. Therefore, the non-local 
branching superprocess is a small particle approximation for the non-local branching 
particle system. Heuristically, ^ gives the law of the migration of the "particles", 0(x, ■) 
describes the amount of offspring born at x G -E by a parent that dies at this point, 
and ipi^r) describes the amount of the offspring born by this parent that are displaced 
randomly into the space according to distributions vr randomly chosen by G{x, dir). Thus 
the locations of non-locally displaced offspring involve two sources of randomness. 

Replacing / in ()3.21|) and ()3.22|) by 6f and differentiating at 6* = we see that the 
first moments of the superprocess are given by 

/ iy{f)Qt{fi,du)=fi{TJ), t>0,feB+iE), (3.23) 
Jm{e) 

where {Tt)t>o is a locally bounded semigroup of kernels on E determined by 

Ttfix) + £ [jj>{y)Tt.sf{y) + m[Tt-sf{y) - m(y,T,_J)]P,(x,d2/)}ds = Ptfm.2A) 
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and m(a;, dy) is the bounded kernel on E defined by 

mix, /) = A(x, /) + / dz/). (3.25) 

Jm{e)° 

In particular, we may define another locally bounded semigroup of kernels {Ut)t>o on E 
by 

Utfix) + ^ {j^m[Ut-sf{y) - m{y, Ut^sf)]Ps{x, dy)}ds = Ptf{x), (3.26) 

which has weak generator G such that 

Gf{x) = Afix) + /3(x)[m(x, /) - fix)], f E ViA), (3.27) 

where A denotes the weak generator of {Pt)t>o- Now the generator B of {Tt)t>o can be 
expressed as 

Bfix) = Gfix) - h{x)f{x), f e ViA). (3.28) 

By a comparison theorem we have Ttf < e"*"*[/t/ for alH > and / G B(E)^. From this 
and ()3.23|) we have 

/ u{f)Qtifi, dz/) < ell''llV(t/t/), t>0,feB+{E). (3.29) 
Jm(e) 

Note that ()3.29j) implies that u \—>- u{l) is a ||6||-excessive function for {Qt)t>o- 

To conclude this section, let us consider briefly the special, and possibly more desir- 
able, case where G{x,d7T) = unit mass at some 7r(x, ■) G Mq{E), that is, the non-locally 
displaced offspring born aX x & E choose their locations independently according to the 
(non-random) distribution Ti{x, ■). In this case, the non-local branching mechanism is 
given by 

^(x, /) = /5(x)[/(x) - C(x, vr(x, /))], xeE,feB+{E), (3.30) 

where 

Cix, z) = d{x)z+ (1 - e"^")ra(x,dM), xeE,z>0, (3.31) 
Jo 

where d G B^{E) and un{x, du) is a bounded kernel from E to (0, oo) with 

/•oo 

m{x) := d{x) + / un{x,du) < 1, x G E. (3.32) 
Jo 

In particular, if ({x, z) = z, we may rewrite ()3.20j) formally as 

jVJix) = AVJix) - 0(x, VJ{x)) + P{x)[n{x, VJ) - VJ{x)], t>0,xeE, 
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with initial condition Vq/ = /. This equation corresponds to a superprocess with underly- 
ing generator A and non-trivial local and non-local branching mechanisms. Alternatively, 
we may also think that the superprocess has underlying generator Af{x) + (3{x)[tt{x, f) — 
f{x)] and only non-trivial local branching mechanism. Since the generator 5 of a general 
Markov process in E is the limit of a sequence of operators of the type P{x)[7i{x, /) — /(x)], 
in principle a superprocess with more general underlying generator A + B and only local 
branching mechanism can be approximated by a sequence of superprocesses with underly- 
ing generator A and non-trivial local and non-local branching mechanisms. Under suitable 
conditions it is also possible to establish convergence of branching particle systems with 
underlying generator A and non-trivial local and non-local branching mechanisms to the 
superprocess with underlying generator A + B and with only non-trivial local branching 
mechanism, which has been done in a particular setting in Gorostiza (1994); see also 
section 6. 



4 Rebirth superprocesses 

We may consider a modification of the branching particle system described in the last 
two sections. Let (^, 7, F) be given as in section 2. A rebirth branching particle system 
with parameters (^,7,F) is described by (2. A), (2.B) and the following 

(4.C) When a particle dies at a point x & E, it gives birth to a random number 
of offspring in E according to the probability kernel F(x,dz/). In addition, the parent 
particle itself is replaced by an extra offspring at site x & E, that is, the parent particle 
is reborn. All the offspring then start to move from their locations. 

Let {Xt : t > 0} be the process defined in the same way as in section 2. Then 
{Xt : t > 0} is still a Markov process with state space N{E). We also have ()2.2|) and 
()2.H|1 . but ()2.4j) is now replaced by 



-ut{x) 



P.{p(0,t)e-^(«')} 



(4.1) 



+ Px 



P(0,s)7(es 



-«t-s(5s)p-i'("t-s) 



N{E) 



F(e.,dz/)]ds}. 



This is equivalent to 

^-utix) ^ P^e-^(«*) -P,|^ 7(^,)e-"'-^(«=)ds} 



or 



I _ Q--^t{x) 



+ Px 
Px{l 

+ p. 



-Ut-s{^s)„-l^(ut~s) 



N{E) 



ds 



,.|^*7(e,)(l-e-"-^(«=))ds} 
(1 - e-"'-=(«=))e-'^("'-^)F(e., du) 



N{E) 



ds 
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+ PJ f Uq [ (l-e-^("'-=))Ffe,dz/)ld4. 

Jo L Jn{E) ^ ^ 

Let Vt{x) = Vt{x, f) be defined by ()2.fij) . Then we have 



VtiX 



7fe) / t;i_.fc)e-^('^*-»)F(e.,dz/) 

N{E) 



ds 



(4.2) 



7(e.) / (l-e-'^("-»))F(e„dz/) 

Af(£;) 



ds 



We now consider a sequence of rebirth branching particle systems {Xt{k) : t > 0} with 
parameters (^, jk, Fk). Define {xl^^^ : ^ > 0} and choose Fk as in section 3 with ak{x) = 
and 7a:(x) = [3k{x). Then ()3.5j) remains vahd if we replace ()3.7p by 



P. {fc(l - e-/(«*)/^)} - P,.{ jju{is)vfl{is)ds] 



+ P. 



{ / k(es) / i;;'Ue.)e^''("-)F,(^„dz/)lds| (4.3) 

Jo L JN{E) ^ 

ft 



+ P. 



4/ \kms) I [l-e-^("-)]F,(e„di/)ld4, 
or equivalently 

vf'\x) + P.[/?fc(es)0fc(es, vf],) + V^fc(e., t;E)]ds = P.A;[1 - e- 
Jo 

where ijjk is given by ()3.10j) and ()3.13j) . and 



(4.4) 



0fc(x, /) = -/(x) / hk{x,'n,l-Ti{f)/k)G{x,dTi). 

'MoiE) 



(4.5) 



Lemma 4.1 /f the conditions of Lemma \3.1\ (i) are fulGUed and if (3k ^ (3 & B~^{E) 
uniformly as k ^ oo, then, for each a > 0, we have (f)k{x,f) — * f{x) uniformly on 
E X B^{E) and the solution v^^\x, f) to l\4.4\) converges boundedly and uniformly on the 
set [0,a] X E X B^{E) of (t,x,f) to the unique bounded positive solution Vtf{x) to the 
evolution equation 



Vtf{x) - / / P{y)ay, Vt.J)Ps{x, dy) ds = PJ{x 



(4.6) 



where (i'r) ^'s defined by \3.15\) and \3.21]) . 
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Based on this lemma, one can show as in section 3 that the finite dimensional distribu- 
tions of {xj:''^ '■ t > 0} under Q^^j converge as A; — oo to those of the process {X^ : t > 0} 
with semigroup {Qt)t>o defined by (j3.22|) and (j4.6|) . Since afc(x) = in the approximating 
sequence, we call {Xt : t > 0} a rebirth superprocess. Note that (j4.6|) is the special form of 
()3.20|) with local branching mechanism (f){x, z) = —(3{x)z, which exactly compensates the 
death factor in the non-local branching mechanism. This observation might be helpful in 
understanding the non-local branching mechanism given by ()3.2H) . 

5 Multitype superprocesses 

In this section, we deduce the existence of a class of multitype superprocesses from 
that of the non-local branching superprocess constructed in section 3 following the ar- 
guments in Li (1992b). Let E and / be two Lusin topological spaces and let ^ = 
{]?, {rjt, at), ^t, P{x,a)} be a right continuous strong Markov process with state space 
E X I. Let ■, ■) and ■, ■) be given by ()3.19p and ()3.31|) . respectively, with x E E 
replaced by {x, a) E E x I. Let l3{-, ■) E E x I and let 7r(x, a, db) be a probability kernel 
from E X I to I. As a special form of the model given in section 3, we have a non-local 
branching superprocess {Xt{dx,da) : t > 0} in M{E x I) with transition probabilities 
determined by 

Q^exp{-Xt(/)} = exp{-^(\4/)}, t > 0, / G B^{E x J), (5.1) 
where Vtf is the unique bounded positive solution to 

Vtf{x,a) + / P(x,a)[(p{Vs, Ois,Vt-sf{Vs, Ois)) + PiVs, Ois)Vt-sf{Vs,Ois)]ds 

Jo 

t 

P{x,a)[l3iVs, as)C{Vs, tts, T^iVs, tts, VtsfiVs, ■)))]ds (5.2) 
= P{x,a)[f{Vt,at)]- 

We may call {Xt : t > 0} a multitype superprocess with type space /. Heuristically, 
{rjt : t > 0} gives the law of migration of the "particles", {at ■ t > 0} represents the 
mutation of their types, 0(x, a, ■) describes the amount of the a-type offspring born when 
an a-type parent dies at x E E, Ci.^i'^i ') describes the amount of the offspring born by 
this parent that change into new types randomly according to the kernel 7r(x, a, d6), and 
(3{x,a) represents the birth rate of the changing-type offspring aX x E E. It is assumed 
that all of the offspring start migrating from the death site of their parent. Note that the 
migration process {rjt : t > 0} and the mutation process {at : t > 0} are not necessarily 
independent. 

Now let us consider a special case which has been studied in the literature. Suppose 
that / = {1, . . . , fc} and for each i E I, r]*-*-* is a right continuous strong Markov process 
in E with semigroup {Pt^^)t>o, belongs to the class given by ()3.19|) and C*^*^ belongs 
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to the class given by (j3.3H) . Let ^ be a right continuous strong Markov process in the 
product space E x I with transition semigroup {Pt)t>Q defined by 

Ptf{x,z) = [ f{y,i)PP{xAy). f e B+{E X /). 
Je 

Let (^((x,^),^) = (x, z). Suppose that G B^(E) and 7c{x,i, ■) is a Markov kernel 
from E X I to I having the decomposition 

k 

i=i 

where p^j\x) > 0, Yl'j=oP^f i^) = 1 6j denotes the unit mass at j G /. Then 
we have a multitype superprocess {Xt : t > 0} in M{E x I) by (jS.lj) and ()5.2|1 . For 
i e I and fi e M{E x I) we define /x^^^ G M{E) by = x {i}). The map 

yU (— > (yU*-^-*, . . . , /i^'^-') is clearly a homeomorphism between M{E x I) and the /c- dimensional 
product space M{E)K Therefore, . . . , X^^^) : t > 0} is a Markov process in the 

space M{E)^, which may be called a k-type superprocess. Clearly, this class of /c-type 
superprocesses coincides with the one defined in Li (1992b). Heuristically, r]^^^ gives the 
law of the migration of the zth type "particles", (f)^^\x,-) describes the amount of the 
ith type offspring born when an ith. type parent dies at point x G E, C^^\x, ■) describes 
the amount of the offspring born by this parent that change into new types randomly 
according to the discrete distribution {Pi\x), . . . ,p^^\x)}, and (5^^\x) represents the birth 
rate of the changing-type offspring at a; G -E. The study of multitype superprocesses was 
initiated by Gorostiza and Lopez-Mimbela (1990); see also Gorostiza and Roelly (1991) 
and Gorostiza et al (1992). 



6 Superprocess-controlled immigration 

By the discussions in the last section, we have a special 2-type superprocess {{x[^\x^^'') : 
t > 0} in M[E)^ with transition probabilities determined by 

exp {-Xi^)(/(^)) - Xf )(/(^))} = exp {-;.(^)(.«) - f^^^\v?^)} , (6.1) 



where vl-^\-) and v^''^-) are defined uniquely by 



,(2) 



vi'\x) + [ {^^'\y,viliy))-vi%iy)) P^'\x,dy)]ds = P,^'^f^'\x), (6.2) 



and 



v?\x) + 



<P^'Hy,vl%iy))P^'\x,dy)]ds = P}'^f'\x). 



(6.3) 
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In particular, if /^^^^ = 0, we have vf^ = and 



(6.4) 



where vl-^\-) is given by 



<P^'\y,vll{y))pW[x,dy)]ds = P}'^f'\x). 



(6.5) 



Thus {XP : t > 0} is a superprocess in M{E) with parameters {r]^^\ (p^^^) . On the other 
hand, by an expression of weighted occupation times, the value in 1)6.111 is equal to 

Q(.(.),.(^)) exp {-X,«(/(i))} exp {-/.(^)(.f ) - X('\vll)ds^ ; 
see e.g. Dawson (1993) and Dynkin (1993, 1994). Then we see that 



Q(M(i),At(2)) 



exp|-X, 



(2) 



X(^) : s > 



exp 



That is, given {Xf ^ : t > 0}, the second coordinate {X^ : t > 0} is a superprocess with 
parameters (r/*^^), 0^^^) and wrf/i immigration controlled by {X^^^ : t >0}. A special class 
of superprocess-controlled immigration processes have been studied in Hong and Li (1999) 
and their relation with multitype superprocesses has been pointed out in Li (2002). 



7 Mass-structured superprocesses 

A multitype superprocess {X((dx, da) : t > 0} with type space / = (0, oo) can be called 
a mass- structured superprocess if we interpret x E E and a > as the coordinates of 
position and mass, respectively. For the mass-structured superprocess, we may consider 
its aggregated process {Yt : t > 0} defined by 

POO 

Yt{dx):= / aXi(da;,da), x E E. (7.1) 
Jo 

Since the integrand on the right hand side is unbounded, {Yt : t > 0} is only well-defined 
under some restrictions. In general, {Yt : t > 0} is not Markovian. Let A be the weak 
generator of {{rit,at) t > 0} and (Tt)t>o the locally bounded semigroup of finite kernels 
on E X I with generator 

Bf{x, a) = Af{x, a) + [3{x, a)[m(x, a)7r(a, /(x, ■)) — /(x, a)] — 6(x, a)f{x, a), (7.2) 

where b{x, a) is the coefficient of the linear term of 0(x, a, z) and m(a;, a) is defined by 
()3.32|) with X E E replaced by {x, a) E E x I. Indeed, ()7.2j) is of the same form as 
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()3.28|) with {x, a) instead of x. By the discussions in section 3, the first moments of the 
sup erpro cess are given by 

Q,{Xt{f)} = fi{TJ), feB^iExI). (7.3) 

In practice, we may have that a newborn offspring is no larger than its parent, which 
corresponds to the assumption that 7r(a, ■) is supported by (0,a). Let H{x,a) = a and 
suppose that H G V^A) is a Ci-excessive function of {{rit,at) : t > 0} for some constant 
Ci > 0. In this case, we have BH{x,a) < (ci + a) and hence H G T>{B) is a 

(ci + ||6||)-excessive function of {Tt)t>Q. It follows from (j7.3j) that 

Q,{X,iH)} < e(^^+II^IIV(/f). (7.4) 

Then we may change the state space slightly and take any cr-finite measure fi on Ex (0, oo) 
satisfying fi{H) < oo as the initial state of {Xf : t > 0}; see e.g. El Karoui and Roelly 
(1991) or Li (1992c). In this case, ^TM implies that Xt{H) < oo a.s. for alH > so that 
()7.1|) defines an aggregated process {Yt : t > 0} with finite measure values. 

A special type of mass-structured superprocess with Markovian aggregated process 
has been studied by Gorostiza (1994). Assume that /3(-,-) = and at = g{t,aQ) for a 
deterministic mapping g{-, ■) from [0, oo) x (0, oo) to (0, oo). Let denote the conditional 
law of {rjt : t > 0} given rjo = x. For / G B^{E), (j5.2|) becomes 

Vtf{x,a)+ P^[(l){7]s,g{s,a),Vt-sf{Vs,g{s,a)))]ds = P^f{r]t). (7.5) 

Since the motion of at = g(t,ao) is deterministic, if Xq is supported by x {a}, then 
Xt is supported by x {g(t,a)} and Yt = g{t,a)Xt. In this case, {Yt : t > 0} is a 
Markov process since the transformation Xt ^ Yt loses no information. For B G ^{E), 
let X^{B) = Xt{B X {g(t,a)}). Then {X^ : t > 0} is an inhomogeneous superprocess 
with cumulant semigroup {V^^)t>r>o defined by V^^f{x) := Vt-rf{x,g{r,a)), which has 
underlying process {rit '■ t > 0} and time- dependent branching mechanism (f){x, g{t,a), ■). 
This gives a representation of the aggregated process in terms of an inhomogeneous su- 
perprocess. A representation of this type was first given by Gorostiza (1994) in the case 
where at = aoe^* for ^ constant c G M. Gorostiza (1994) obtained the process as high 
density limit of a sequence of branching particle systems where the mass of each offspring 
is equal to that of its parent multiplied by a fixed positive constant factor, and the mass of 
any particle does not change during its lifetime, realizing in a particular case the program 
mentioned at the end of section 3. 

8 Multilevel superprocesses 

Multilevel superprocesses arise as limits of multilevel branching particle systems. In a two 
level system, objects at the higher level consist of non-trivial sub-populations of objects 
at the lower level and both lower level and higher level objects can branch. A lower level 
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object consisting of a population can be described by a measure on some space S. We 
can then view a two level system as a multitype system with / = M{S)°, the space of 
non-trivial finite Borel measures on S. Non-local branching is natural in this context. For 
example, at the particle level the offspring of a second order object consisting of a set of 
particles could consist of a subset of the particles or include more than one copy of the 
original particles. 

To make this precise, we may let 5" be a topological Lusin space and {at : t > 0} be the 
Markov process with state space M{S)° obtained by killing a superprocess at its extinction 
time. Then {Xt : t > 0} is a Markov process with state space M{E x M{S)°), which 
can be called a m,v,ltilevel superprocess generalizing the model of Dawson and Hochberg 
(1991), Dawson ct al (1990) and Wu (1994). For the multilevel process, it is also natural 
to study the aggregated process {Yt : i > 0} defined by 

Yt{AxB):^[ [ a{B)Xt{dx,da), A e B{E), B e B{S). (8.1) 

J A J M{S) 

To illustrate the possibihties of non-local branching, consider the case in which E is 
a singleton. In this case, we may view {Xt : i > 0} as a superprocess with state space 
M{M{S)°). A possible non-local branching mechanisms is obtained by taking 7r(/x, dv) to 
be the law of 



N 



J2^z^^ (8-2) 



where A?^ > 1 is an integer-valued random variable, and {Zi,Z2,---} are i.i.d. random 
variables in S with distribution fi{S)~^ fi{-) . That is, the offspring of a level two object 
H e M{S) is a single point measure with the same total mass as /i and its location is 
selected randomly according to the empirical measure of a sample from the normalized 
parent distribution. 

Another possibility is given by 

7r(//,di/) = 5^,(di/), (8.3) 

where B e B{S) and //s e M{S) is defined by //b(A) = //(A n B). In this case the 
offspring is a level two object in which only level one individuals falling in the set B (Z S 

are present. 

In the case in which ii^ is a countable set, we may interpret the M{E x M(5')°)-valued 
process {Xt : i > 0} as a population in a sequence of islands. The £^-coordinate tells in 
which island the M(S') °-valued objects {at : t > 0} in the first level are located. The 
non-local branching is given by n{x,iJ,,dp) = 7r{/j,du), which only acts on the M(S)°- 
coordinate at the higher level. Jumps in E described by {rjt '■ t > 0} correspond to 
the independent migration of "clans" (families of the lower level) between the islands. 
Suggestively, we may call {Xt : t > 0} a stepping stone type superprocess. 

Properties and applications of multilevel superprocesses involving local branching have 
been studied extensively in the literature; see e.g. Dawson and Hochberg (1991), Dawson 
et al (1990, 1994, 1995), Etheridge (1993), Gorostiza (1996), Gorostiza et al (1995), 
Hochberg (1995), Wu(1994) and the references therein. 
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9 Age-reproduction-structured superprocesses 



Let ii^ be a Lusin topological space and let ^ = {f2, {rit,at,Ot), ^ , J^t, P{y,a,z),'y} be a 
Borel Markov process with state space E x M"*" x N^, where 7 is a terminal time; see 
Sharpe (1988, p. 65). We assume that both at and 6t are non- decreasing processes. Let 
/3(-, •, ■) G B{E X M+ X N+) and let C(-, ■) be given by (nOT|l with x e E replaced by 
(x, a,z) G -E X M+ X As a special form of the models given in sections 4 and 5, we 
have a rebirth multitype superprocess {Xt : t > 0} in M{E x x N^) with transition 
probabilities determined by 

Q^exp{-Xi(/)} = exp{-/i(V,/)}, t > 0, / G x M+ x N+), (9.1) 

where Vtf is the unique bounded positive solution to 

Vtf{y,a,z) - / P(y,a,z)[P{Vs,as,Os)C{Vs,as,Os,Vt_J{r]„0,0))l{^^^^}]ds (9.2) 
Jo 

= Piy,a,z)[fiVt,O:t,0t)l{a,<j}]- 

It is not hard to check that the first moments of the superprocess are given by 

Q,{Xtif)} = KTtf), f e B+{E X M+ X N+), (9.3) 
where (Tt)t>o is a semigroup of bounded linear operators on B^{E x M+ x N+) defined by 

Ttf{y,a,z) - / P(y,a,z)[P{Vs,(ys,0s)m{r]s,a„9s)Tt-sfiVs,0,0)l{a^^<-y}]ds (9.4) 
Jo 

= P{y,a,z)[f{Vt,at,Ot)l{a,<'y}], 

where m(-, ■, ■) is given by ()3.32|) with x E E replaced by (x, a, z) G x M"*" x N^. Using 
{Tt)t>o we may rewrite ()9.2|1 into 

Vtf{y,a,z)+ [ Ts[(3mVt-sf-a-r,;Vt^sf{;0,0my,a,z)ds = Ttf{y,a,z). (9.5) 

In the case at = ao + t, we may call {Xt : t > 0} an age-reproduction- structured 
superprocess. Heuristically, rjt represents the location of a "particle", at its age and 9t 
the number of its offspring born in the time interval {t — at,t]. At each branching time, 
the particle gives birth to a random number of offspring whose motions start from the 
branching site and whose ages and reproduction numbers start from zero. The particle 
does not disappear at its branching times, it is removed from the population only when its 
age exceeds the lifetime 7. An interesting limit theorem for age-reproduction-structured 
branching particle systems was proved in Bose and Kaj (2000) which leads to the super- 
process in the special case where E' is a singleton and rjt = r]o. (Compare ()9.5|) and their 
equation (2.8).) 
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